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on : 

We present a new timing solution for PSR B1534+12, based on coherently-dedispersed ob- 
servations at Jodrell Bank and, recently, Arecibo. The new data have resulted in improved 
measurements of the post-Keplerian timing parameters, including the orbital period deriva- 
tive, Pi,. At present, the poorly-known distance to the pulsar limits the precision of the 
measurement of the intrinsic Pb, and hence the strength of the test of general relativity that 
results from this binary system. By assuming that general relativity is the correct theory of 
' gravity, we may invert the test and find an improved value of the pulsar distance. 

Oh; 

o : 

, 1 Introduction 

Double-neutron-star binary pulsars in close, highly eccentric orbits have long provided the best 
. ^ | strong- field tests of the predictions of gravitational theories. The timing analysis of pulsar 

^ ■ signals permits the measurement of five Keplerian orbital elements as well as a number of 

post-Keplerian (PK) orbital parameters. The PK parameters can be determined usinga theory- 
independent procedure in which the masses of the two stars are the only unknownsu Each of 
the PK parameters depends on the masses in a different way; consequently, if any two of them 
are measured, the relevant parameters of the two-body system are fully determined within any 
gravitational theory. If three or more PK parameters can be measured, a test of the gravitational 
theory results from the overdetermined system. 

For the first binary pulsar, PSR B1913+16, the PK parameters Co (rate of advance of perias- 
tron), 7 (time dilation and gravitational redshift) and Pb (orbital period derivative) have been 
measured and found to be in excellent agreement with the predictions of general relativityBi'Bi 
PSR B1534+12, discovered in 1990^ is a comparable system, with an eccentric 10.1-hour orbit. 
PSR 131534+12 is significantly brighter than PSR B1913+16, and its pulse profile has a narrow 
peak, permitting timing measurements of nearly five times greater accuracy. Because the orbit is 
nearly edge-on as viewed from the Earth, the PK parameters r and s (the "range" and "shape" 
of the Shapiro time delay) are measurable, in addition to w, 7 and Pb- The resulting overde- 
termination of the orbit leads to a non-radiative test of gravitation theory in the strong-field 
regime, complementing the Cj-^f-Pb test for PSR B1913+160 

Previously published timing measurements of PSR B1534+12 were made with the Arecibo 
305 m telescope through early 1994pl when the telescope went out of normal service for a major 



upgrading; and from 1994 through 1997 with the 43 m telescope at Green Bank, West Virginia 
and the 76 m Lovell telescope at Jodrell Bank, England! With the recent reopening of Arecibo, 
we have been able to conduct a new series of observations of this pulsar in 1998 July, and improve 
on our previously reported timing results; in particular, we have reduced the uncertainty on the 
observed value of the PK parameter Pf, by 35%, to approximately 5% of the expected GR value. 



2 Observations 



A complete discussion of the observing procedures and instrumentation may be found in ref. 9. 
We only note here that over the years, new instrumentation was developed to remove more 
perfectly the pulse-smearing effects of interstellar dispersion. While the 1990-1994 Arecibo ob- 
servations and the Green Bank observations used filterbanks to divide the observing bandpass 
into small channels, the later observations were carried out with an improved "coherent dedis- 
persion" system, which sampled the raw telescope voltages and then removed the effects of 
dispersion using a software filter. The more recent observations therefore yield significantly 
better timing precision. 

All observations involved folding the pulse signal over several minutes at the predicted 
topocentric pulse period to produce a total-intensity profile. Cross-correlation with a stan- 
dard template, using a least-squares method in the Fourier transform domain, yielded a time 
offset This offset was added to the time of the first sample of a period near the middle of 
the integration, thereby yielding an effective pulse arrival time. A different standard template 
was used for each observing system and frequency; they were made by averaging the available 
profiles over several hours or more. Uncertainties in the TOAs were estimated from the least 
squares procedure, and also from the observed scatter of the TOAs within 30 minutes of each 
one. Each observatory's local time standard was corrected retroactively to the UTC timescale, 
using data from the Global Positioning System (GPS) satellites. 



3 The Timing Model 

A pulse received on Earth at topocentric time t is emitted at a time in the comoving pulsar 
frame given by 



T = t - t + A c - D/f + A RQ + A EQ -A 
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Here to is a reference epoch and Ac is the offset between the observatory master clock and the 
reference standard of terrestrial time. The dispersive delay is D/f 2 , where D = DM/2.41 x 10~ 4 , 
with dispersion measure DM in cm~ 3 pc, radio frequency / in MHz, and the delay in seconds. 
Finally, Arq, Aeq, and A$q are propagation delays and relativistic time adjustments for effects 
within the solar system, and Ar, Ae and As are similar terms accounting for phenomena within 
the pulsar's orbitHa'O The orbital A terms are defined by: 
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These are written in terms of the eccentric anomaly u and true anomaly A e (u), and the time 
dependence of u, which are related by: 
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A e (u) = 2arctan 
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At a given time t, then, the propagation delay across the pulsar orbit is calculated by a model 
which incorporates ten parameters implicitly defined in the above equations: five Keplerian 
parameters (x, u, Tq, P&, e) and five PK parameters (lj, P^ 7, r, s). These quantities, in 
conjunction with a simple time polynomial to model the spin of the pulsar and with astrometric 
parameters to model the propagation of the signal across the solar system, constitute the free 
parameters to be fit in the theory-independent timing model. 

In a particular theory of gravity, the five PK parameters can be written as functions of the 
pulsar and companion star masses, mi and mo. and the well-determined Keplerian parameters. 
In general relativity, the equations areH'E 
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Here the masses mi, 1712, and M = mi + m 2 are expressed in solar units, and we use the 
additional shorthand notations s = sini and T Q = GM Q /c 3 = 4.925490947 /js, where i is the 
angle between the orbital angular momentum and the line of sight, G the Newtonian constant 
of gravity, and c the speed of light. 



4 Arrival Time Analysis 

We used the standard tempo analysis softwareil together with the JPL DE200 solar-system 
ephemerist3 to fit the measured pulse arrival times to the timing model with a least-squares 
technique. Arbitrary offsets were allowed between the different data sets to allow for frequency- 
dependent changes in pulse shape and slight differences in the standard profile alignments. 
To partially compensate for small systematic errors, uncertainties in TOAs were increased (in 
quadrature) by 2.9, 17.1, 20.5, 7.0, and 6.0 fis, respectively, when calculating weights for the 
1990-94 Arecibo, Jodrell Bank, Green Bank data, 1998 Arecibo 430 MHz and 1998 Arecibo 
1420 MHz sets. 

Results for the astrometric, spin, and dispersion parameters of PSR B1534+12 are presented 
in Table |]. We fit the data to two models of the pulsar orbit. The theory-independent "DD" 
modelEJ treats all five PK parameters defined in §H| as free parameters in the fit. Alternatively, 
the "DDGR" mo assumes general relativity to be correct and uses equations || through |l2| 

to link the PK parameters to M = mi + m 2 and m 2 ; consequently it requires only two post- 
Keplerian free parameters. 

Table |2| presents our adopted orbital parameters. Uncertainties given in the table are approx- 
imately twice the formal "1 <r" errors from the fit; we believe them to be conservative estimates 
of the true 68%-confidence uncertainties, including both random and systematic effects. The 
Keplerian orbital parameters include the period P^, projected semi-major axis x = a\smi/c, 



Table 1: Astrometric, spin, and dispersion parameters for PSR B1534+12. Figures in parentheses are uncertainties 
in the last digits quoted, and italic numbers represent derived quantities. 



Right ascension, a (J2000) 


15 h 37 m 09^959952(16) 


Declination, 5 (J2000) 


11° 55' 55"6561(3) 


Proper motion in R.A., u n (masyr -1 ) 


1.5(2) 


Proper motion in Dec, (masyr" 1 ) 


-25.6(3) 


Parallax, tt (mas) 


< 1.5 


Pulse period, P (ms) 


37.9044404878550(14) 


Period derivative, P (10 -1 ^) 


2.42250(3) 


Epoch (MJD) 


48778.0 


Dispersion measure, DM (cm _3 pc) 


11.6152(16) 


DM derivative (cm _3 pcyr _1 ) 


-0.00004(53) 


Galactic longitude I (deg) 


20.0 


Galactic latitude b (deg) 


47.8 


Composite proper motion, [i (masyr -1 ) 


25.6(3) 


Galactic position angle of // (deg) 


238. 7(4) 



eccentricity e, longitude of periastron u, and time of periastron To. These quantities are followed 
by the measured post-Keplerian parameters relevant to each of the two models. The DDGR 
solution includes a parameter called "excess -Pf,," which accounts for an otherwise unmodeled 
acceleration resulting from galactic kinematics. 

The best estimates of the masses of the pulsar and its companion come from the DDGR 
solution. We find the masses to be mi = 1.344±0.002 M and m 2 = 1.335±0.002 M . For the 
sake of comparison Table § lists, in italic numbers, computed PK parameter values corresponding 
to the measured masses in the DDGR fit. The fitted and derived parameter values are in accord, 
indicating good agreement of the theory-independent solution with general relativity. 

Figure [l] shows the post-fit residuals for the 1990-1994 Arecibo 1400 MHz data, the Jo- 
drell Bank data, and the recent Arecibo 430 MHz data, plotted as functions of date. Even 
within a single data set, the TOA uncertainties can vary by factors of three or more because of 
scintillation-induced intensity variations. We have not attempted to show these differences in 
data quality by means of error bars in the residual plots. 



5 Discussion 



Before the observed P^ can be compared to the theoretical value, we must apply a correction 
which accounts for the relative acceleration of the center-of-mass frame of the binary pulsar 
system and the solar system bary center. An expression for the most significant bias, which 
arises from galactic kinematic effects, is derived in ref. 3. This bias can be written as the sum 
of terms arising from acceleration toward the plane of the Galaxy, acceleration within the plane 
of the Galaxy, and an apparent acceleration due to the proper motion of the binary system: 
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Here a z is the vertical component of galactic acceleration, I and b the galactic coordinates of 
the pulsar, Rq and vq the Sun's galactocentric distance and galactic circular velocity, [i and d 



Table 2: Orbital parameters of PSR B1534+12 in the DD and DDGR models. Figures in parentheses are 
uncertainties in the last digits quoted, and italic numbers represent derived quantities. 





DD model 


DDGR model 


Orbital period, Pb (d) 


0.42073729933(3) 


0.42073729933(2) 


Projected semi-maior axis, x (s) 


3.729464(3) 


3.7294638(5) 


Eccentricity, e 


0.2736775(5) 


0.2736774(2) 


Longitude of periastron, u (deg) 


267.44738(10) 


267.44744(9) 


Epoch of periastron, Tn (MJD) 


48777.82595097(6) 


48777.82595096(6) 


Advance of periastron, uj (degyr -1 ) 


1.755794(19) 


1.75580 


Gravitational redshift, 7 (ms) 


2.071(6) 


2.067 


Orbital period derivative, (P b ) obs (10~ 12 ) 


-0.131(9) 


-0.1924 


Shape of Shapiro delay, s 


0.983(8) 


0.9762 


Range of Shapiro delay, r (/js) 


6.3(1.3) 


6.6 


Total mass, M = mi + 7772 (M & ) 




2.67845(4) 


Companion mass, 777,2 (M&) 




1.344(2) 


Excess P b (10~ 12 ) 




0.061(9) 



the proper motion and distance of the pulsar; we use the short-hand notation (3 = d/ Rq — cos I. 
The pulsar distance can be estimated from the dispersion measure, together with a smoothed- 
out model of the free electron distribution in the Galaxylll This model yields d ~ 0.7 kpc for 
PSR B1534+12, with an uncertainty of perhaps 0.2 kpc. At this distance we estimate a z /c = 
(1.60 ± 0.13) x 10~ 19 s _1 El Following ref. 3, we assume a flat galactic rotation curve and take 
vo = v\ = 222 ± 20kms~ 1 and Rq = 7.7 ± 0.7 kpc. Then, summing the terms in equation (|i~3|) 
and multiplying by Pb, we find the total kinematic correction to be 

(iV) gal = (0.038 ± 0.012) x 10" 12 . (14) 

The uncertainty in this correction is dominated by the uncertainty in distance, which is only 
roughly estimated by the dispersion-measure model. 

Our measurement of the intrinsic rate of orbital period decay is therefore 

/ . \ obs / . \ gal n 

(PfeJ - (AJ = (-0.169 ± 0.015) x 10~ 12 . (15) 

Under general relativity, the orbital period decay due to gravitational radiation damping, (Pb) , 
can be predicted from the masses m\ and 777-2 (eq. |i~0| ), which in turn can be deduced from the 
high precision measurements of uj and 7. The expected value is 

GR, 

Pb) = -0.192 x 10~ 12 . (16) 

Although the measured value in the pulsar center-of-mass frame differs from this prediction 
by 1.5 standard deviations, it can be brought into good agreement by increasing the pulsar 
distance to slightly over 1 kpc. Stated another way, we can assume that GR is the correct 
theory of gravity, measure the "excess -FV' for the system as described above and presented in 
Table [2L and then invert equation 13 to determine the pulsar distancet3 Figure |2| shows the 



relation between the two quantities; using this approach we obtain d = 1.08 ± 0.15 kpc ^0070 
confidence limit). The uncertainty is dominated by the measurement uncertainty of (Pb) obs , 
rather than uncertainties in the galactic rotation parameters or the acceleration a z . Continued 
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Figure 1: Post-fit residuals versus date for (a) Arecibo 1400 MHz, (b) Jodrell Bank and (c) Arecibo 430 MHz 

data. 



timing of this system should lead to a much more precise distance measurement; in fact, the 
recent Arecibo observations have already reduced the uncertainty by 35% over our previously 
published value! 

We note that the timing solution provides a second, independent constraint on the distance. 
The upper limit on parallax, ir < 1.5 mas (Table constrains the distance to d > 0.67 kpc. 
While the parallax distance has less precision than the kinematic distance, it is reassuring that 
these measurements are in agreement. 

An accurate distance for PSR B1534+12 is of considerable interest because this system, 
along with PSR B1913+16, is of prime importance in estimating the rate of coalescence of 
binary neutron-star pairs in a typical galaxy. Early estimates of the inspiral rate used much 
smaller distances for PSR B1534+12, including 0.4 kpc@ 0.5 kpc@ and 0.7 kpcpH leading 
to a low estimate of its intrinsic luminosity. Our distance, compared to that estimated from 
the dispersion measure, increases the volume per B1534+12-type object by a factor of 2.5 to 4, 
depending on the unknown scale height of such systems in the Galaxy; this lias helped lead to 
improved estimates of the neutron-star-neutron-star coalescence rate&EiE3£3 



5. 1 Test of Relativity 

This pulsar provides the second test of general relativity based on the to, 7, and parameters 
of a binary pulsar system, and is the first double-neutron-star binary to permit significant 
measurements of the Shapiro-delay parameters r and s. The left-hand sides of equations (|| O) 
represent measured quantities, as specified for this experiment in the "DD" column of Table |2[ 
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Figure 2: Relation between measured "excess P;," and pulsar distance for the PSR B1534+12 system. Dashed 
lines indicate the 68% confidence ranges of the measured values. 

If GR is consistent with the measurements and there are no significant unmodeled effects, we 
should expect the five curves corresponding to equations (|8|-|T2|) to intersect at a single point 
in the m\-m2 plane. A graphical summary of the situation for PSR B1534+12 is presented in 
Figure ||, in which a pair of lines delimit the 68% confidence limit for each PK parameter (a 
single line for Cb, whose uncertainty is too small to show). A filled circle at mi = 1.335 Mq, 
m2 = 1.344 Mq marks the DDGR solution of Table [2L and its location on the Co line agrees well 
with the measured DD values of 7 and s. These three parameters therefore provide a unique 
test of the purely quasi-static regime of general relativity. We have already noted that the DD 
value of Pb is slightly too small when corrected to the dispersion-estimated distance. However, 
as discussed above, this discrepancy can be removed by invoking a larger distance to the pulsar. 
Finally, the value of r, although presently little better than a 20% measurement, also agrees 
with its expected value. This pulsar thus provides a convincing second demonstration of the 
existence of gravitational radiation; furthermore, all its timing parameters are in good accord 
with general relativity theory. 
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